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Abstract 



The coherent fi-e conversion rates in various nuclei are calculated for general lepton 
flavor violating interactions. We solve the Dirac equations numerically for the initial state 
muon and the final state electron in the Coulomb force, and perform the overlap integrals 
between the wave functions and the nucleon densities. The results indicate that the 
conversion branching ratio increases for a light nucleus up to the atomic number Z ~ 30, 
is largest for Z = 30 - 60, and becomes smaller for a heavy nucleus with Z > 60. We also 
discuss the uncertainty from the input proton and neutron densities. The atomic number 
dependence of the conversion ratio calculated here is useful to distinguish theoretical 
models with lepton flavor violation. 
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1 Introduction 



The observation of lepton flavor violation (LFV) is one of the most interesting signals 
beyond the Standard Model (SM). The charged-lepton LFV processes such as the [i — > e-f 
decay and the \x-e conversion in muonic atoms can occur in many promising candidates 
beyond the SM, although the simple seesaw neutrino model does not induce experimen- 
tally observable rate for the \x — > ey process. For example, sleptons in the supersymmetric 
(SUSY) extension of the SM and bulk neutrinos in the higher dimensional models generate 
LFV processes through one- loop diagrams P], 0]. In the R-parity violating SUSY models, 
additional LFV interactions exist at the tree level The branching ratios of the LFV 
processes have been calculated in many models in the literature, especially for supersym- 
metric grand unified theories (SUSY-GUTs) [f|, H and a SUSY model with right-handed 
neutrinos 0. It was shown that fi — > ej and fi-e conversion branching ratios can be close 
to the experimental bounds in these models. 

There are on-going and planned experiments for the \x — > ej and fi-e conversion 
searches. For the \x — > cy branching ratio, the present upper bound is 1.2 x 10~ n from 
the MEGA collaboration 0. A new experiment is under construction at PSI aiming 
for a sensitivity of 10~ 14 ||. For the /i-e conversion, an upper bound for the conversion 
branching ratio is 6.1 x 10 -13 || reported by the SINDRUM II experiment at PSI. Now 
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SINDRUM II is running with gold (Au) targets. The MECO experiment at BNL (TC 
planned in order to search for fiAl — > eAl with a sensitivity below 10~ 16 . In future, further 
improvements by one or two orders of magnitude are considered for both fi — > &y and fi-e 
conversion processes in the PRISM project ]TT] at the new 50 GeV proton synchrotron 
constructed as a part of the JAERI/KEK joint project. 

In order to compare the sensitivity to the LFV interaction in various nuclei, precise 
calculation of the fi-e conversion rate is necessary. There have been several calculations 
of the conversion rate. Weinberg and Feinberg calculated in the case that the conversion 
occurs through the photonic interactions (/i-e-7 vertex) ||12|| . In the calculation, they used 
several approximations in which the muon wave function was taken to be constant in a 
nucleus and the outgoing electron was treated as a plane wave. The plane wave treatment 
of the electron is a good approximation only for light nuclei because the effect of Coulomb 
distortion on the electron wave function is large for heavy targets. The non-photonic in- 
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teraction case was studied by Marciano and Sanda ]T3|]. Shanker improved the calculation 
by solving the Dirac equations for the muon and electron wave functions in the electric po- 



tential of a nucleus [ 14 1 . The calculation was carried out for all the interactions including 
the photonic and four-fermion operators in the effective Lagrangian, but the treatment of 
the photonic dipole operator was incomplete because he used an approximation that off- 
shell photon exchange was replaced by the four-fermion interaction. Recently, Czarnecki 
et al. presented calculation in which the off-shell photon is properly treated as an electric 
field in a nucleus and listed the values of the conversion rate for aluminum (Al), titanium 
(Ti), and lead (Pb) targets [IE . The transition rate to the ground state of a nucleus as 
well as excited states are calculated in Refs.|TB|, |T7| ]. 

In this paper, we evaluated the \x-e conversion rates for nuclei of a wide range of 
atomic numbers by the method of Czarnecki et al. We took into account all the operators 
for the fx-e transition. For any types of operators, the results of our calculation indicate 
a tendency that the conversion branching ratio is larger for the nuclei with moderate 
atomic numbers than that for light or heavy nuclei. Although the tendency is the same, 
there are significant differences in Z dependences of the conversion rate for various LFV 
couplings. The experiments in various nuclei are therefore useful for model discrimination 
because each theoretical model predicts different Z dependences. The conversion rate 
depends on the input proton and neutron densities for each nucleus. Although the proton 
density is well measured by electron scattering, there is uncertainty in the determination 
of neutron densities. We estimate the uncertainty from these input parameters. Based on 
neutron density distribution determined from proton scattering experiments performed in 
1970's and pionic atom experiments, the conversion rate changes by 20% - 30% for heavy 
nuclei. The ambiguity is shown to be significantly reduced for lead (Pb) by the improved 
determination of the neutron density from a new proton scattering experiment. 

This paper is organized as follows. In section ^ we present a formula of the conversion 
rate with the most general effective Lagrangian for LFV. The results of our calculation 
and the estimation of the uncertainty are shown in section || In section f|, we summarize 
this paper. The model parameters of the nucleon density functions in nuclei and the muon 
capture rate in nuclei are listed in Appendices |A| and [B|, respectively. 
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2 Formula of /x-e conversion rate 



In this section, we present a method of the conversion rate calculation. We solve the 
Dirac equations for the muon and electron in the initial and final states, respectively, and 
obtain transition amplitudes by integrating the overlap of both wave functions. 

We start with the most general LFV interaction Lagrangian which contributes to the 
fi-e transition in nuclei 



where Gp and are the Fermi constant and the muon mass, respectively, and A^r 
and g's are all dimensionless coupling constants for the corresponding operators. Our 



conventions are = d»A„ - d^, a*" = (i/2)[j», 7"], P L = (l-j 5 )/2, P R = (l + 7s)/2, 



and the covariant derivative is defined as = d^—iQeA^, where Qe (e > 0) is the electric 
charge (Q = — 1 for the electron and the muon). The size of each coupling constant 
depends on the interaction of the new physics in which the lepton flavor conservation is 
violated. There are two types of amplitudes for photonic transition (/x-e-7), namely the 
monopole and dipole fi-e transitions. In the above expression, the monopole transition is 
converted to the vector-vector interaction assuming that the momentum dependences of 
the form factors are negligible. 

The initial state in the /i-e conversion process is the Is state of the muonic atom, and 
the final electron state is the eigenstate with an energy of m M — e^, where et, is the binding 
energy of the Is muonic atom. Both wave functions in the initial and final states can be 
determined by solving the Dirac equations in the electric field of the nucleus. The Dirac 
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equation in the central force system is given by |T 
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where W and V(r) are the energy and potential, respectively, m 8 is the reduced mass of 
the electron or the muon, cr are the Pauli matrices, and the orbital angular momentum I is 
defined by I = — ir x V. Since the operator K and the ^-component of the total angular 
momentum j z commute with the Hamiltonian, two eigenvalues of these operators, — k 
and (jl, represent the quantum numbers of the wave functions describing this system. We 
denote the wave function as follows: 



0(r)x{J(M) 



(4) 



where Xk is the normalized eigenfunction of (cr ■ I + 1) and j z such as 



1 i-2-k 

dcosO 
i Jo 



xTx" K . = <V<W . (5) 



The total angular momentum j is related to k as k = ±(j + 1/2). With the notation of 
m (r) = rc/(r) and U2(r) = rf(r), the Dirac equation for the radial function is given by 



d ( u\ 
dr \ u 2 



—k/t W — V + rrii 

— (W — V — m-i) k/t 



"2 



(6) 



The initial muon state corresponds to the quantum numbers of fi = ±1/2 and K — — 1 
with a normalization of 



d 3 *^* (x)#(x) = l. 



(7) 



The final electron state is one of the states in the continuum spectrum. Our normalization 
convention is taken as 



d 3 ^^*(x)^,^(x) = 27r5^5 KK ,5(W - W) 



(8) 
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The conversion rate uj conv in this normalization is simply written by the square of the 
fi-e conversion amplitude M, taking the spin average of the initial muon and summing 
over the final states of the electron. From the effective Lagrangian ([I]), M is obtained as 
follows: 

4G F 



M 



V2 

a 

V2 



J d?x (m,A*J^a"Pp R ^ + m,Al^%a^P L ^) (N'\F aP \N) 
9LS {q )^%PR^ + Qrs^wPl^) (N'\qq\N) 



3 x 

q=u,d,s 

A e ) td „/,(m) i „ , ,„7>( e ) jd /at'If;, 



+ (^l^^^l^M^Sfrf) (N'\gj a q\N) 
+ (gLA^SrPL^ +9RA {q) €§l a P R 4' L s) (N'\qj a J 5 q\N) 



\ (9LT (q) € ( ^ aP P R ^ + 9RT {q) €%a a "P L ^) (N'\qa aP q\N) 



(9) 



where (N'\ and \N) are the final and initial states of the nucleus, respectively. 

Hereafter, we concentrate on the coherent conversion processes in which the final 
state of the nucleus is the same as the initial one. The fraction of the coherent process 
is generally larger than incoherent one approximately by a factor of the mass number of 
the target nuclei. In this case, the matrix elements of (N\qj 5 q\N) , (N \ q^alsq \N), and 
(N\qa a pq\N) vanish identically, and (N\qq\N) and {N\q^ a q\ N) can be expressed by the 
proton and neutron densities (p^ and p< n >) in nuclei as follows: 

(N\qq\N) = ZGf p) p (p) + (A - Z)Gf n) p {n) , (10) 

2Zp^ + (A - Z)pV> for q = u, 

(N\q 1 °q\N}={ Zp^ + 2(A-Z)p^ for q = d, (11) 

for q = s, 

(N\qfq\N) = (i = 1,2,3). (12) 

Here we introduce coefficients G^' pS> and G^' n ^ for scalar operators. These are evaluated 
to be G^' P) = Gf n) = 5.1, Gf p) = Gf n) = 4.3, and G^ p) = G^ n) = 2.5 in Kosmas et 
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al. [19|. We assume that the proton and neutron densities are spherically symmetric and 
normalized as follows: 



drinr 2 p^ n \r) = 1 . 



(13) 



The final formula of the conversion rate can be written as follows: 



2Gi 
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(14) 



A* R D + g®Sto + ~9^S^ + g®vV> + ~gtlv {n) 
AID + g®SV + ~ 9 P S S^ + fifW y W + § Wy(n) 

Al and Ar are given in Eq.(|Ij). The coupling constants g's in Eq.(|H|) are defined as 

-gilns = E^'^b), (15) 
9ls,rs = J2 G s' n) 9ls,rs( q ) , (16) 



~(p) 

9LV,RV = ^9LV,RV(u) + gLV,RV(d) 



~(n) .) 
9LV,RV ~ 9LV,RV(u) + ^9LV,R(d) • 

We also introduced the following overlap integrals: 
4 f°° 

D = Jo drr ^~ E ^ 9 ^^ + f* 9 > 



1 Z* 00 



1 f 00 

2^2 L drr2Zp[P) ^^ + ^^ 



n ) > 



(17) 
(18) 

(19) 
(20) 
(21) 
(22) 
(23) 



where the functions g e etc. are defined in the Is muon wave function and k = ±1 electron 
wave functions as follows: 

±1/2/ a ,\ 



1 if;{r)xT'\e,< 



(24) 



€:^ 2(e, M» = ( i^t'l''^ I ■ (26) 



In the above expression, we have neglected the electron mass, so that g+ and /+ are 
related to g~ and /~ as g+ = if~ and i/+ = By integrating the Maxwell equations, 
the electric field E(r) is determined as 

E(r) = ^ / rVV')^ • (27) 



o 



Once the proton and neutron densities are given, one can calculate the electric field 
E(r) by Eq. (p7|) and the electric potential V(r) by 



^(r) = -e / E(r')dr' . (28) 



The wave functions are then obtained by the Dirac equation Eq.(^|), and the /i-e conversion 
rate is calculated by Eq.flU]). 

3 Numerical results 

In order to evaluate Eq.(|Hp, we need proton and neutron densities and the muon and 
electron wave functions. We first discuss proton and neutron densities and show the 
wave functions of muon and electron. Then, we present numerical results of the overlap 
integrals Eqs.(^) - (|23|) and the conversion rates for various nuclei. 

3.1 Distribution of protons and neutrons in the nucleus 

We have used proton densities determined from electron scattering experiments. In the 
past, the charge density distribution of a nucleus was analyzed assuming some func- 
tional form of the proton distribution such as the two-parameter Fermi model, the three- 
parameter Fermi model, and the three-parameter Gaussian model. More recently, with 
improvement of experimental data, model-independent methods are used to extract de- 
tailed information on the density distribution. Examples are the Fourier-Bessel expansion 
and the sum of Gaussian functions. The proton density is very precisely determined when 
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the model- independent analysis is carried out. We use the charge density listed in Ref. |2C] . 



We adopt the results of the model-independent analysis when the data are available. More 
on charge distribution is described in Appendix [A|. 

The determination of the neutron distribution is not as easy as that of the proton 



distribution [21]. There are several ways to determine the neutron density in the different 
region. Pionic atoms provide a probe of the neutron density on the periphery of the 
nucleus. In a certain level of the pionic atom the pion is absorbed by the nucleus. Since 
the strong interaction between the pion and the nucleus changes the energy and the width 
of this level, we can obtain information on the neutron density in the nucleus from the 
analysis of the atomic X-ray spectrum. Scattering experiments on the nucleus by strong 
interacting particles such as the proton, the alpha particle, and the charged pion are also 
useful to determine the neutron density. Recently, an experiment with antiprotonic atoms 
was carried out to determine the proton and neutron density distribution in the periphery 



region of various nuclei [22 



In this paper we use the following three methods to evaluate the fi-e conversion ratio: 

1. First, we take the proton density from electron scattering experiments given in 
Appendix |A] and assume the neutron density is the same as the proton density. 
For light nuclei this is a good approximation because the number of neutrons and 
protons are approximately equal and the conversion rates do not depend on the 
details of the neutron distribution, (method |TJ) 

2. Second, we employ the neutron distribution obtained from the pionic atom. We use 



the results of the analysis of Ref. [E3[, where the proton and neutron distributions 



are given in terms of the two-parameter Fermi function, (method |2|) 

3. Finally, we use the neutron distribution obtained from the polarized proton scatter- 
ing experiment. The analysis was carried out for carbon (C) |[24| , titanium (Ti) f25| , 



nickel (Ni) [23], zirconium (Zr) an d lead (Pb) where proton and neutron 
density are given in the literature. We also estimate the uncertainty due to the error 
of the neutron distribution from the scattering experiment based on Refs. [^, |27|, [28| . 
(method |3|) 

The first method gives precise evaluation for the overlap integrals D, and V^. On 
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the other hand, the neutron density is necessary in order to determine and V^ n \ 
Since both methods |2| and |3| use strong interacting particles as probes to the neutron 
density, the calculation suffers from ambiguities associated with the pion/proton-nucleon 
interaction. The method [| provides information on the size of the nucleus in a wide range 
of atomic numbers. On the other hand, we can determine in the method |3] the profile of 
the neutron distribution inside the nucleus for selected nuclei. We calculate the conversion 
rate in two methods, and first discuss ambiguities in each method and then compare the 
results. 

3.2 Numerical evaluation of the overlap integrals 

In this subsection, we first show an example of the muon and electron wave functions 
obtained by solving the Dirac equation Eq.(fJ), and present the result of numerical calcu- 
lation of overlap integrals D, S^' n ' and V^ ,r ^ defined by Eqs.(19) - (p3|). 



3.2.1 Wave functions of the initial and final states 

The muon and electron wave functions are evaluated by solving the Dirac equation (|2]) 
with the electric potential given by Eq.(^Sj). Ignoring the recoil of the nucleus which is 
of the order of m^/M/v, where is the nucleus mass, one finds that the energy of the 
out-going electron in Eq.(FJ) is equal to the muon mass minus the binding energy. As 
an example, we show the muon and electron wave functions in titanium (Ti) nucleus in 
Figs.|H and [2|. We can see in Fig.jl] that the lower component m 2 in the muon wave function 



is much smaller than the upper component U\. However, as pointed out in Ref. [[14]], its 
effect on the conversion rate is sizable for heavy nuclei. The overlap integrals are evaluated 
using these wave functions. 

3.2.2 Method 1 

We present in Table [I] the results under the assumption p n = p p , where p p is taken from 
Ref. |2(J . We show in Fig.^j the Z dependence of the integrals. We omitted the points for 



gg 6 Er, ylpTa, and yg 7 Au from this figure since these data are either obtained from quite 



old experiments or not well-established GO] . We see that the overlap integrals increase as 



functions of Z for light nuclei up to Z ~ 30, and saturate or decrease for heavy nuclei. 
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Table 1: The overlap integrals in the unit of m 5 J 2 are listed. The proton distribution in 
the nuclei are taken from Ref. (see also Appendix [A]), and neutron distribution are 
assumed to be same as that of protons (method [I] in subsection |3~T| ). 
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Table 1: (Continued). 
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Figure 1: The normalized wave function of a muon in the titanium (Ti) nucleus is plotted. 
The solid line and dashed line represent Ui(= rg~) and u 2 (= rf~) components, respec- 
tively. The horizontal axis is the distance between the nucleus and the muon in the unit 
of l/m M . The unit for the wave function is taken to be JrfL 



Although the tendency is the same, each overlap integral has different Z dependence, 
especially for heavy nuclei. For example, the scalar (S^ p ' n ') and the vector (V^ p,n - ) ) type 
integrals are almost equal in light nuclei (Z < 30), whereas the vector like integral is 
larger by a factor of 1.5 - 2 than the scalar one for heavy nuclei. This difference is due 
to the the relativistic effects of the muon wave functions which are significant in heavy 
nuclei. In fact, the scalar and vector overlap integrals in Eqs. (|20|) and ( p2|) [Eqs. (|21|) and 
(p3])1 are exactly the same, if we ignore the small component of the wave function /~ For 
D we can show that D/(8e) ~ ~ is satisfied for light nuclei. 

The qualitative feature of the Z dependence of the overlap integrals can be understood 
from the following consideration. When we adopt an approximation by Weinberg and 



Feinberg [12|, where the muon wave function is replaced by the average value and the 
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Figure 2: The normalized wave function of an electron in titanium (Ti) nucleus is plotted. 
The solid line and dashed line represent Ui(= rg~) and u 2 {= rf~) components, respec- 
tively. The horizontal axis is the distance between the nucleus and the electron in the 
unit of l/m M . 



electron wave function is treated as a plane wave, the formula for the overlap integrals in 
Eqs.([20|) and ( p2|) are given by 
1 



y(p) ^ g(p) 



57T 



p ■ 



Here F p is the form factor defined by 



dr4irr 2 p (p) — 



and ((pfj) is the average value of the muon wave function in the nucleus given by 



POO 

( ) 2 = / drAnr 2 (gl + fl) p 
Jo 



{p) _AnvyZin 



(29) 



(30) 



(31) 



In the last expression, we have introduced Z e ^ which is the effective charge for the muon 
in the Is state. We show Z e s in Fig.f|. Since the muon wave function enter into the inside 
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Figure 3: The overlap integrals are plotted as functions of the atomic number Z. Neutron 
density distributions are assumed to be equal to charge density distributions (method |I] 
in subsection |3~T| ). 



of the nucleus, Z e s thus does not increase linearly with respect to Z. In Fig.[|, we show 
the form factor F p calculated based on method [l|. There we see that F p is a decreasing 
function and suppressed for heavy nuclei. These two properties explain the Z dependence 
of the overlap integrals. 

3.2.3 Method 2 

We present in Table and Fig.^| the values of the overlap integrals in the method 0, 
in which the input nucleon distributions are obtained from the analysis of the pionic 
atom experiments |[2^]. In this method, the proton and neutron density are assumed to 
be a two-parameter Fermi model defined by Eq. (fH|) in Appendix |X[ Using the energy 



shift and the decay width of the pionic atom, the neutron density is determined together 
with parameters in the pion-nucleon optical potential. In Ref.^3], it is assumed that the 
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Figure 4: The effective charge for the muon in the Is state is plotted as a function of the 
atomic number Z. 



diffuseness of the neutron density is the same as that of the proton density, z n = z p , so 
that the output parameter is the radius of the neutron density. For our calculation, we 
use the neutron matter parameter i? n [mean] in Table 4 of Ref. |23|| . The Z dependence 
shown in Fig|| is similar to that seen in Fig.[3]. 

In order to estimate the ambiguity within this method, we have calculated and 
y( n ) by changing the radius of the neutron density c„ within the error given in the Ref. |23| . 



We show the results for 13 Al, j^pCe, and gffPb corresponding to light, intermediate, and 
heavy nuclei in Table ^. The error of c n includes the statistical error as well as the 
systematic error estimated from the outputs in two different parameterizations of the 
optical potential |§. We see that and change ±2% for j 7 3 M, ±(6 - 8)% for 
If Ce, and ±(7 - 11)% for |> 8 Pb. 
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Nucleus 


D 






S [n) 


yW 


9 r 


0.0166 


0.0071 


0.0072 


0.0089 


0.0090 


f?Na 


0.0260 


0.0111 


0.0114 


0.0128 


0.0131 


ffMg 


0.0299 


0.0128 


0.0132 


0.0126 


0.0131 


?s 7 Al 


0.0357 


0.0153 


0.0159 


0.0163 


0.0169 


£Si 


0.0421 


0.0181 


0.0188 


0.0173 


0.0180 


32 q 
16° 


0.0529 


0.0227 


0.0238 


0.0221 


0.0233 


f 8 Ar 


0.0628 


0.0268 


0.0284 


0.0310 


0.0330 


40 p„ 
20^ a 


0.0778 


0.0333 


0.0355 


0.0319 


0.0341 


26 re 


0.110 


0.0464 


0.0508 


0.0503 


0.0555 


63 p 
29 ^U 


0.124 


0.0521 


0.0579 


0.0585 


0.0654 


gGe 


0.138 


0.0576 


0.0651 


0.0704 


0.0802 


lAs 


0.141 


0.0585 


0.0665 


0.0690 


0.0792 


^ 2 Ce 


0.188 


0.0698 


0.0934 


0.0844 


0.117 


ifAu 


0.167 


0.0523 


0.0859 


0.0610 


0.108 


i 8 Pb 


0.162 


0.0495 


0.0838 


0.0575 


0.107 


i 9 Bi 


0.163 


0.0495 


0.0846 


0.0579 


0.108 



Table 2: Same as Table [I], but here the results of the analysis of the pionic atom 
experiment are used for the distribution of the neutrons p3f (method [| in subsection 
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Figure 5: The form factor F p is plotted as a function of the atomic number Z. 



3.2.4 Method 3 

The neutron densities for selected nuclei were determined from the 800 MeV polarized pro- 
ton elastic scattering experiments performed at LAMPF in the late 1970's. In Refs. [[24|, |25|] , 
proton and neutron densities are given assuming the three-parameter Fermi or Gaussian 
model. Errors of the neutron distributions are estimated in a model independent fash- 
ion in Refs. p6| |27|j . More recently, the determination of the neutron density have been 
improved for lead (Pb) based on a new experiment and a new analysis |28|. [29| . 

We calculate the overlap integrals using neutron distribution given in the above ref- 
erences. Table ^ and Fig.[7] show the results for the five nuclei from the Refs. p3|, |2"5fl . 
In order to evaluate uncertainty from the neutron distribution determined by the proton 
experiments, we take several examples where the error of the neutron distribution is ex- 
plicitly given in the literature. We calculated and within the uncertainty of the 
neutron distribution. For example, we present the error band of the neutron distribution 
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Figure 6: Same as Fig.^|, but here the results of the analysis of the pionic atom experiment 
are used for the distribution of the neutrons [|33| (method |2] in subsection 



in the nucleus g2 8 Pb in Fig.|8] |[28|| . Similar error bands based on older experiments are 



given in Ref. [26] for foCa, ||Ni, and 5Q 6 Sn, and in Ref. [27| for 82 8 Pb. The overlap in- 
tegrals evaluated using the minimum and maximum values of the envelope are shown in 
Tables | and |. Since the proton distribution is not explicitly given except the case of 
Pb in Ref. [p8fl , we use the proton distribution listed in Table [7] for other cases. The 



208 
82 

errors of and amount to a few percent for the light nuclei such as 20 Ca and 
2gNi. For I^Pb, we can see drastic improvements in the determination of and V^ n ' 



based on the new measurement. According to Ref.||28||, the reduction of the errors in the 
neutron distribution is mainly due to higher statistical accuracy in the new experiment. 
This implies that improvement in the proton scattering is important to determine neutron 
densities more precisely and reduce the ambiguity in the calculation of the /x-e conversion 
rate for heavy nuclei. 
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Nucleus 


c n [fm] 






£A1 


3.09 ± 0.08 
6.00 ± 0.10 
6.86 ± 0.09 


0.0163 =F 0.003 
0.0844 t 0.0067 
0.0575 T 0.0066 


0.0169 =F 0.003 
0.117 =F 0.008 
0.107 =f 0.008 



Table 3: Error of the overlap integrals associated with the input value of the neutron radii 
determined by the pionic atom experiment. The larger (smaller) values of S^ 1 ' and 
correspond to smaller (larger) value of c n . 



Nucleus 


D 






g(n) 




6 ^ 

iTi 

INi 

90 7 r 

i 8 Pb 


0.0074 

0.0870 

0.130 

0.176 

0.156 


0.0032 
0.0371 
0.0548 
0.0715 
0.0457 


0.0032 
0.0399 
0.0605 
0.0844 
0.0812 


0.0031 
0.0462 
0.0606 
0.0841 
0.0712 


0.0032 

0.0495 

0.0667 

0.100 

0.122 



Table 4: Same as Table [I], but here the results of the analysis of the proton scattering 
experiments are used for the neutron density distribution [|4], (method || in subsection 

nu 



3.2.5 Comparison of the results in three methods 

Comparing Tables [I], @, and |], one finds that the overlap integrals of the light nuclei agree 
with one another within a few percent. For heavy nuclei, the method [l] is not necessarily 
a good approximation for S^ n ' and V^ n \ However, we see that the results are consistent 
with the values in the method |3] within 10%. Comparing the method |^ and ^ the pionic 
atom method gives generally smaller values by 10 — 20% than the analysis based on the 
proton scattering for the intermediate and heavy nuclei. 

We should note that the present analysis with pionic atoms assumes the two-parameter 
Fermi model with z p = z n for neutron distributions. In Ref . p2 , the neutron distribution 
is analyzed from an antiprotonic atom experiment, namely the nuclear spectroscopy anal- 
ysis of the antiproton annihilation residues and the measurements of strong-interaction 
effects on antiprotonic X-rays. The authors compared two neutron distributions of neu- 
tron skin type (z n = z v and c n > c p ) and halo type (z n > z p and c n = c p ), and concluded 
that the halo type distribution was favored. 

In order to illustrate how the results depend on the assumption of the neutron distri- 
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Figure 7: Same as Fig|| but here the results of the analysis of the proton scattering 
experiments are used for the neutron density distribution (method ^ in subsection 

HI 



bution, we calculate the and in the halo type distribution with fixed values of 
the mean square radius. For lead (Pb), (V^) increases by 22% (12%) compared to 
the skin type distribution, so that the results are very close to the values in the method ||. 
This indicates that the overlap integrals are sensitive to the neutron distribution inside 
the nucleus, not only global quantities such as the mean square radius. In the method 
|2], the pionic atom data provide us with information on the neutron distribution mostly 
in the peripheral region, and therefore the interior neutron distribution is given by an 
extrapolation based on the assumed functional form. Since the skin type distribution is 
assumed in Ref.|23], there may be errors associated with this assumption. 



We have shown that in the method |] the ambiguity of and is reduced to 
a few percent even for heavy nuclei such as lead (Pb) if we use the results of the new 
experiment. The neutron distribution derived by experiments depend on how we treat the 
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Figure 8: The error envelope of the neutron density distribution for Pb nucleus [f28 . 



scattering of a proton and a nucleus, so that there might be uncertainties in the evaluation 
of and associated with validity of the scattering theory. The improvement of 



the scattering theory would be necessary to give more realistic error estimation |30 
3.3 Numerical evaluation of conversion rate 

It is now straightforward to evaluate the conversion rate through Eq. (0) once a theoretical 
model with LFV is specified. In this subsection, we present a fi-e conversion branching 
ratio for various types of LFV interactions in order to show the possibility to discriminate 
the different models through the Z dependence. We also compare our results with the 
existing calculations for the case where the photonic dipole operators are non-vanishing. 
We consider the following three cases: 



1. The photonic dipole operators Al and/or Ar are non- vanishing. The fi-e conversion 
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S (n) 


T T ( T) \ 

y( n i 


Ref. 


2gCa Minimum 
Maximum 


0.0331 
0.0338 


0.0352 
0.0359 


m 


2gNi Minimum 
Maximum 


0.0584 
0.0592 


0.0644 
0.0651 


m 


5Q 6 Sn Minimum 
Maximum 


0.0958 
0.104 


0.120 
0.128 


[26 


g2 8 Pb Minimum 
Maximum 


0.0679 
0.0789 


0.120 
0.131 






Table 5: Maximum and minimum values of the overlap integrals for the neutron density 
distribution, which is changed within the error envelope. 





D 


S (p) 




S (n) 




Ref. 


g2 8 Pb Minimum 
Maximum 


0.163 
0.163 


0.0493 
0.0493 


0.0845 
0.0845 


0.0675 
0.0697 


0.0119 
0.0121 


[28 



Table 6: Maximum and minimum values of the overlap integrals for the neutron density 
distribution for g2 8 Pb based on the analysis in Ref. 



branching ratio is given by 

_ 2G 2 F D\\A L \ 2 + \A R \ 2 ) 



B 



OJ, 



UJ. 



capt 



(32) 



capt 



where w C a,pt is the muon capture rate. For convenience, we list the capture rate in 



Appendix B RTJ 



2. The scalar operators gRs(d) and/or gis(d) are non- vanishing. The /x-e conversion 
branching ratio in this case is given by 



B 



fj,N-*eN 



2GUGf p) S^ + Gf n) S^n\g LS{d) \ 2 + \g RS{ 4' 

^capt 



(33) 



3. The vector operators guv(u) and gLV(u) satisfies gnv(u) = —2gRV(d) 7^ and/or 
9lv(u) — —2gLV(d) 7^ . The />e conversion branching ratio in this case is given by 

2G 2 v^(\~gil\ 2 + \^l\ 2 ) 



B 



(34) 



capt 
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Figure 9: The /i-e conversion ratios for the typical theoretical models are plotted as func- 
tions of the atomic number Z . The solid, the long dashed, and the dashed lines represent 
the cases that the photonic dipole, scalar, and vector operator dominates, respectively. 
Proton and neutron distribution are taken according to method [1] in subsection |3J], and 
the conversion ratios are normalized by the conversion ratio in aluminum nuclei (Z = 13). 

The first case appears as a good approximation in SUSY models for many cases, especially 
in SO(10) SUSY GUT models [|J and in SUSY models with right-handed neutrinos ||. 
The second case is realized in some cases of SUSY models with R-parity violation ||. The 
third case corresponds to the situation where the monopole form factors give dominant 
contributions in the /i-e-7 transition. The /i-e conversion branching ratio are shown for 
three cases in Figs|5| (method [TJ), [LT] (method |2|), and [TTJ (method ^). In these figures the 
branching ratios are normalized by the value for aluminum evaluated by the method [T], 
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Figure 10: The fi-e conversion ratios for the typical theoretical models evaluated by 
the method |^ in subsection |TT]. The branching ratio is normalized by B^ N ^ eN (Z = 13) 
evaluated by the method |l|. 

namely 

Dipole: B, lN ^ eN {Z = 13) = 9.9 (\A L \ 2 + |^| 2 ) , (35) 

Scalar: B^ eN (Z = 13) = 1.7 x 10 2 {\g L s(d)\ 2 + \g R s(d)\ 2 ) , (36) 

Vector: B^ eN {Z = 13) = 2.0 (\~g^\ 2 + | 2 ) . (37) 

We can see that, for all three types, the branching ratio increases as Z for Z < 30, are 
largest for 30 < Z < 60, and decreases for Z > 60. It is also seen that the conversion 
ratios have large differences in heavy nuclei depending on the three types of interaction. 
From this property we may be able to distinguish models beyond the SM through several 
experiments with different targets. 

In order to see improvements of the present method from older calculations, we com- 
pare three different approximations for the case where the photonic dipole operators are 
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Figure 11: The fi-e conversion ratios for the typical theoretical models evaluated by 
the method |^ in subsection |TT]. The branching ratio is normalized by B^ N ^ eN (Z = 13) 
evaluated by the method |l|. 

non-vanishing: namely, our calculation, Weinberg-Feinberg approximation, and the ap- 
proximation by Shanker. For this purpose we define the ratio R(Z) = B^N-^eN / 'B(p — > 
ej), where the \i — > branching ratio £>(/i — > ej) is given by 3847t 2 (|Al| 2 + |A^| 2 ). The 
present method thereby gives 
r 2 n 2 

R W = iqt*, • (38) 
1927r^u; capt 

In the Weinberg-Feinberg calculation, the relativistic effects and the Coulomb distortion 
were ignored [12|]. We define the conversion branching ratio Bj^^ eN and the ratio of 
ratios R WF (Z) in the Weinberg-Feinberg approximation by the following formula: 

B^ eN = 8 -^* 3 Z* s ZF;(\A L f + \A R f)^, (39) 

7T ^capt 

R {Z) ~ - e 7 ) " 48vr^ capt ' (40) 

26 



Notice that these are not exactly the same as the formula given in the original paper 
because they used approximate formulae for the capture rate and the form factors for 
the general photonic transition. Shanker improved the Weinberg- Feinberg formula taking 
into account the relativistic effects and the Coulomb distortion. In his approximation, 
the branching ratio and the ratio of ratios for the dipole photonic interaction are given 
by 

B^ eN = l28e 2 GlV^(\A L f + \A R f)—, (41) 

^capt 



R s(z) = rvv^v c F ,-- (42) 



B s ^ eN _ e 2 GlV^ 
B(jji -> erf) 37r 2 u; ca pt 

We present our R(Z), R WF (Z), and R S {Z) in FigJT2. Here we used the proton density 



in the method [I] and the muon capture rate uj cap t from the experiments pi |. We see that 
the three quantities have similar Z dependence: they range from 0.002 to 0.006, and are 
largest for Z = 30 - 60. The values of R WF (Z) and R S (Z) are larger than our R(Z) by 
30% for Z > 50. We have reproduced with a good accuracy the result by Czarnecki et 
ai, where they evaluated R[Z) for aluminum (Al), titanium (Ti), and lead (Pb) nuclei. 
Kosmas obtained in Ref.]l7|] the result that R(Z) is a monotonically increasing function, 
but he did not take into account the Coulomb distortion effect. We could indeed obtain 
the increasing function of R(Z) by ignoring this effect. Thus the Coulomb distortion 
effect of the wave function is important in the calculation of the conversion rate for heavy 



nuclei, as noted by Shanker [14]. 

4 Summary 

We have calculated the coherent fi-e conversion rate for general LFV interactions for 
various nuclei. We have used updated nuclear data for proton and neutron distributions 
and taken into account the ambiguity associated with neutron distribution. We gave 
the list of the overlap integrals in Tables p], |2|, and ^ for various nuclei, from which 
we can calculate conversion rates for general interactions for LFV using Eq . ([14]) . We 
also investigate the Z dependence of the conversion rate. We see that the branching 
ratio increases for the light nuclei such as Z < 30, are largest for Z = 30 - 60, and 
decreases for heavy nuclei with Z > 60. Although this tendency of Z dependence is 
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Figure 12: The /i-e conversion branching ratio divided by the \x — > decay branching 
ratio for method [I] is plotted as a function of atomic number Z . The solid line (R(Z)), the 
long-dashed line (_R WF (Z)), and the dashed line (R S (Z)) represent the results of our cal- 
culation, the Weinberg-Feinberg formula, and the Shanker's approximation, respectively. 

the same for different types of coupling constants, there are significant differences in 
Z dependence of branching ratios. We show that the ambiguity in the calculation of 
the overlap integrals associated with proton densities (D, and V^) is quite small 
because the charge densities of nuclei are well known. On the other hand, the overlapping 
integrals and contain uncertainly from the neutron distribution, especially for 
heavy nuclei. We have estimated and with several inputs. Using the neutron 
density distribution determined from proton scattering experiments performed in 1970's 
and pionic atom experiments, we showed that the conversion rate changes by 20% - 30% 
for heavy nuclei. Adopting the improved neutron density distribution determined by the 
new proton scattering experiment, we found that the ambiguity is significantly reduced 
down to a few percent. Because the main ambiguity for heavy nuclei is associated with 
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the neutron density, it will be possible to make a precise prediction if we can determine 
the neutron density with improved analysis and experiments. 

The results of our calculation are useful to choose the appropriate target nuclei for 
future experiments for the /i-e conversion search. In addition, it may be possible to 
identify the theoretical models beyond the Standard Model through the Z dependence of 
different interactions when the signal of the \i-e conversion is experimentally observed. 
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A Proton and neutron densities in nuclei 



We introduce models of nucleon densities in nuclei and list the values of parameters of 
these models used in the calculation. 



We used one of the following models for each nucleus [20 



1. Harmonic oscillator model (HO): 



Pp(n) ( 



Po 



1 + a 



cxp 



T \ 2 

a 



(43) 



2. Two-parameter Fermi model (2pF): 



Po 



r - Cp(n)) I 



(44) 



3. Three-parameter Fermi model (3pF): 

p (l + w p{n) r 2 1 c 2 p(n) ) 



P P (n){r) 



1 + exp[(r - c p{n ))/z p{r 



(45) 
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4. Two-parameter Gaussian model (2pG): 

Po 



Pp(n) ( 



l + exp[(r 2 -c 2 )/ 



p(n) > I Z p(n) 



(46) 



5. Three-parameter Gaussian model (3pG): 

( v p (l + w p(n) r 2 /c; (n) ) 
PpW(rj "l + exp[(r 2 -cJ (n) )/4 n) ] 



(47) 



Here c p ( n ), 2 p ( n ), and u> p ( n ) are the model parameters and po is the normalization constant. 
We also used the following model-independent analysis for several nuclei. 



6. The Fourier- Bessel expansion (FB): 

yj avjoivnr/R) for r<R 

V 

for r > R 



Pp(n) ( 



(48) 



where a v are the coefficients, R is the cutoff radius, and the function j (z) = sin z/z 
is the spherical Bessel function of the zeroth order. 



7. The Sum of Gaussian expansion (SOG): 



Pp(n) ( 



exp 



r-Rj 

7 



+ exp 



r + Rj 

7 



(49) 



where 



Ai 



ZeQi 



2tt 3 / 2 7 3 (1 + 2i? 2 /7 2 ) 



(50) 



We list the model and its parameters used in calculation in Table 0. We do not list 
parameters for FB and SOG there; see Ref . [f20fl . 

B Muon capture rate in nuclei 



We list in Table § the muon capture rates uj capt which are used in our calculation 
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Nucleus 


Model 


c p or a p (fm) 


z p (fm) or a 


w p 




SOG 


— 


— 





3 Li 


HO 


1.77 


0.327 







HO 


1.791 


0.611 







HO 


1.69 


0.811 





fc 


FB 


— 


— 







3pF 


2.570 


0.5052 


-0.180 


8 6 


FB 













2pF 


2.59 


0.564 


_ 




3pF 


2.791 


0.698 


-0.168 


1Z o 


3pF 


2.791 


0.698 


-0.168 


? 3 7 A1 


FB 










3pF 


3.340 


0.580 


-0.233 


31p 
i ^-^ 

10 


FB 











32 o 
16 


FB 











ffci 


3pF 


3.476 


0.599 


-0.10 




FB 


— 


— 







3pF 


3.743 


0.585 


-0.201 


2QCa 


FB 











48 T ; 

22 


FB 











51y 

no V 


2pF 


3.91 


0.532 





iCr 


FB 











5fMn 

zo 


2pF 


3.89 


0.567 





26 


3pG 


3.505 


2.325 


0.380 


iCo 


2pF 


4.158 


0.575 


— 


iNi 


FB 








iCu 


2pF 


4.218 


0.596 




SSZn 


3pG 


3.664 


2.425 


0.342 


iGe 


FB 








80c 

3 4 oe 


2pG 


4.622 


2.365 




88 C r 
38 M 


FB 








89v 
39 1 


2pF 


4.86 


0.542 




90 7-, 


FB 









Table 7: The model parameters of the proton density functions are listed. These val- 
ues are extracted from Refs.|2(J. The abbreviations HO, 2pF, 3pF, 2pG, 3pG, FB, and 
SOG represent the harmonic oscillator model, the two-parameter Fermi model, the three- 
parameter Fermi model, the two-parameter Gaussian model, the three-parameter Gaus- 
sian model, the Fourier-Bessel expansion, and sum of Gaussian, respectively. We do not 
list here parameters for FB and SOG; see Ref.[^C] . 
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Nucleus 


Model 


c p or a p (fm) 


z p or a (fm) 


w p 


l?Nb 


2pF 


4.953 


0.541 







FB 


— 








jj°pd 


FB 


— 










2pF 


5.314 


0.571 





115T 


2pF 


5.357 


0.563 







3pG 


5.110 


2.619 


0.292 




2pF 


5.32 


0.57 


— 


i 8 Ba 


3pG 


5.3376 


2.6776 


0.3749 


139t 


2pF 


5.71 


0.535 






2pF 


5.839 


0.569 




f 8 Sm 


FB 








6 4 Gd 


FB 








165 TT 

67 n ° 


2pF 


6.12 


0.57 




166 PY 
68 ^ r 


3pF 


5.98 


0.446 


0.19 


Tel 


2pF 


6.38 


0.64 




7 fw 


2pF 


6.51 


0.535 






2pF 


6.38 


0.535 




28*Hg 


FB 








205 Tl 
81 11 


SOG 








l 8 Pb 


FB 








i 9 Bi 


FB 








232 Th 

90 111 


2pF 


6.851 


0.518 




238 tt 
92 u 


2pF 


6.874 


0.556 





Table 7: (Continued). 
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Nucleus 


^capt (,1U S J 


Nucleus 




4 He 


0.000336 




5.569 


7 3 Li 


0.0018 


75 A 
33 As 


6.104 


iBe 


0.0074 


IfSr 


7.02 


11 T, 

5 B 


0.0219 


39 1 


7.89 




0.0388 


QOrv 


8.66 




0.0693 


l?Nb 


10.36 


8 u 


0.1026 


Q9n r 


9.614 


1 Q-r-i 

9 r 


0.229 


1in-n i 

iiU Pn 
46 ru 


10.00 


?oNe 


0.200 


Is Cd 


10.61 


Q Q TV T 

ffNa 


0.3772 


li0 Tn 
49 111 


11.40 




0.4841 


120 c 
50 1311 


10.44 


97 A i 

13 Al 


0.7054 


191 rn 


10.21 


28 Q ' 
14 D1 


0.8712 


i 8 Ba 


9.94 


15 r 


1.1185 


57 ^ a 


10.71 


'39 n 

16° 


1.352 


58°Ce 


11.60 


35 m 


1 OOO 


142 KM 
60 Wd 


Iz.oU 




1.30 


f 8 Sm 


12.22 


19^ 


1.849 


64 Gel 


11.82 


40 p„ 
20^ a 


2.557 


67 nu 


12.95 


48 Ti 

22 11 


2.59 


166 T7< 

68 ^ 


13.04 


51y 
23 v 


3.069 


Tel 


12.86 


52 p r 
24 ^ r 


3.472 


7 fw 


12.36 


|Mn 


3.857 


7 fAu 


13.07 


iFe 


4.411 


205T1 
81 11 


13.90 


|?Co 


4.940 




13.45 


iNi 


5.932 


i 9 Bi 


13.10 


63 p 
29 


5.676 


9?Th 


13.1 


30 ^ 


5.834 


238 tt 
92 u 


12.4 



Table 8: The total capture rates used in calculation are listed |3lj] 
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